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0^. Abstract 

(~| ' The main background for the supersymmetric stop direct production search comes from Stan- 

dard Model ti events. For the single-lepton search channel, we introduce a few kinematic variables 
to further suppress this background by focusing on its dileptonic and semileptonic topologies. All 
^ , are defined to have end points in the background, but not signal distributions. They can substan- 

ce ' tially improve the stop signal significance and mass reach when combined with traditional kinematic 

variables such as the total missing transverse energy. Among them, our variable Mj^ has the best 
overall performance because it uses all available kinematic information, including the on-shell mass 
of both W^'s. We see 20%-30% improvement on the discovery significance and estimate that the 8 
^L? ' TeV LHC run with 20 fb^^ of data would be able to reach an exclusion limit of 650-700 GeV for 

direct stop production, as long as the stop decays dominantly to the top quark and a light stable 
neutralino. Most of the mass range required for the supersymmetric solution of the naturalness 
problem in the standard scenario can be covered. 



1 Introduction 

A main goal of the Large Hadron Collider (LHC) experiments is to understand electroweak symmetry 
breaking. In the Standard Model (SM), it is achieved by the vacuum expectation value (VEV) of 
a scalar Higgs field. However, a fundamental scalar field receives quadratically divergent radiative 
contribution to its mass-squared and suffers from the hierarchy problem. One of the most promising 
solutions to the hierarchy problem is supersymmetry (SUSY) which introduces a superpartner to every 
SM field, so that the quadratically divergent corrections to the Higgs mass-squared can be canceled 
between SM particles and their superpartners. Supersymmetry has been extensively searched for at 
colliders, and so far we have not found any evidence for it. The latest LHC search results constrain the 
masses of the gluino and (light generation) squarks in the minimal supergravity [U El [3] or constrained 
minimal supersymmetric standard model (MSSM) [3] to be greater than about 1 TeV OE]. At face 
value, it may imply a serious fine-tuning of the electroweak scale if SUSY exists. However, as the 
largest radiative correction to the Higgs mass in the SM comes from the top quark loop, only the top 
superpartners (stops) need to be light enough to cancel the top loop contribution [318]. The gluino 
and first two generation squarks can be heavier than 1 TeV without a naturalness problem, at least 
at one-loop level [3 EJ El [IHl HH 12] • Therefore, searching for the top superpartners at the LHC offers 
the most important test of whether SUSY provides a natural solution to the hierarchy problem. 

Many third generation squark searches at the LHC rely on gluino production, with subsequent 
decay to stops or sbottoms [HI [HI [151 [161 ITT] . This is because the production cross section for gluino 
is much larger than for direct stop or sbottom production, as long as the gluino mass is not much 
heavier. However, since the naturalness of the electroweak breaking scale does not require gluino to 
be light enough to be copiously produced at the 7 or 8 TeV LHC, a more robust stop search would 
only rely on direct stop pair production. Li this paper we focus on the stop search in this channel in 
a standard i?-parity conserving SUSY scenario. Here the stop decays to a top quark and the lightest 
supersymmetric particle (LSP), which is assumed to be a neutralino. Although a light neutralino 
is not required by naturalness, it avoids the stable charged particle problem and provides a natural 
candidate for dark matter. The signal we are looking for is tt + EJ^^^^ where the missing transverse 
energy ii^™'^^ comes from the pair of neutralino LSP's which escape the detector. We also assume that 
the mass difference between the stop and the LSP is substantially larger than the top quark mass. 
Otherwise the signal will strongly overlap with the SM backgrounds, which would need some different 
search strategies [TSl [T9l [20] [211 [22] . 

In fact, the ti + EJ^^^^ signal also occurs in many other extensions of the SM that include a 



dark matter candidate. It is quite natural to mitigate the hierarchy problem with a relatively light 
top "partner" that decays to the top quark and the dark matter particle. This is possible if the 
partner is also charged under the symmetry that protects the stability of the dark matter particle. 
Examples are the little Higgs models with T-parity | !231 [Ml [25l I26j . models with the exotic fourth 
generation and dark matter [27], models with gauged baryon and lepton numbers [28], and so on. 
Consequently, the same search applies to many different models, but the mass-reach depends on 
each model's top partner production cross section. Studies of the ti + i?™^*^*^ signal for new physics 
have been performed by many groups in various (fully hadronic, single- lepton) channels in recent 

years [SSlETlEniEolEIlEaEallMlES]. 

The ATLAS collaboration at the LHC has done such a search in the single-lepton channel based 
on 1.04 fb~^ of data [36]. In the single-lepton channel, one W from the top decays leptonically and 
the other W decays hadronically. Requiring one lepton in the final state suppresses QCD multijet 
backgrounds tremendously while still retaining a significant W branching fraction. The final state 
signal consists of four (or more) jets (including two 6-jets), one lepton, and missing transverse energy. 
Besides the standard transverse momentum px and pseudo-rapidity |r/| requirements for each object, 
the two main variables used for separating the signal and backgrounds are the missing transverse 
energy EJ^^^^ and the transverse mass Mt constructed from the lepton and EJ^^^^ |36j . The signal 
events are expected to have large E^^^^ from the top-partner decays, along with a neutrino from the 
leptonic W. A hard cut on Ej^^^^ very effectively reduces the SM backgrounds. The cut on Mt removes 
the backgrounds where the E^^^^ is mostly due to a single neutrino from a W decay because the Mt 
distribution has an end point at My\r in such cases. The existing ATLAS study focused on a fermionic 
top partner and can exclude this partner's mass up to 420 GeV [36j . There was no sensitivity to the 
SUSY stop with this limited amount of data because of the much smaller cross section for the scalar 
particles. Last year's run already delivered more than 5fb~^ of data. This year, the LHC is expected 
to deliver even more luminosity at a higher center of mass energy of 8 TeV. Given the importance of 
the stop (and other top-partner) search, it is desirable to extend the mass reach using current and 
future data. 

The single-lepton channel analysis of the ATLAS top partner search paper |36| found that the 
largest background remaining after the their cuts on E^'^^^ and Mt is the dileptonic ti. In these 
background events, both W^s decay leptonically, but one of the leptons is not reconstructed, is outside 
the detector acceptance, or is a r lepton (which may be misidentified as a jet). Each event contains 
at least two neutrinos that can produce a large E™^^^ and also make it easier to pass the Mt cut. The 
additional jets come from QCD initial state radiation (ISR). The next-to-largest background comes 



from the semileptonic ti and T^+jets. The other backgrounds are small after the E^^^^ and Mt cuts. 
To improve the search reach, we designed kinematic variables to identify tt backgrounds, focusing on 
its decay topology. We find that the signal significance and mass reach can indeed be substantially 
improved with the help of these variables. This paper is organized as follows. In the next section we 
discuss several such variables. In section [3l we compare the performances of the basic set of cuts and 
cuts including the new variables, identifying an economical set. 

2 Kinematic Variables for the ti Backgrounds 

We study the LHC search for the pair-production of stops, pp — )• tit J, with t^ — )• t + xiO We focus on 
the signal's one-lepton decay channel, in which one top quark decays leptonically t — )■ W^b — )• l^vnb 
and the other one decays hadronically t — t- bjj (also the other way around). The signal contains four 
jets, one lepton, and missing transverse energy. According to the latest ATLAS ti+ E^^^^ search [36], 
the largest SM background after the E^^^^ and Mj- cuts is ti in the dileptonic channel with one lost 
lepton and two additional jets from ISR that fake the hadronic W. In this section we try to identify 
some kinematic variables, based on these background event topologies. 

Before we discuss the new kinematic variables, we first examine the distributions of the signal 
and main backgrounds in some traditional kinematic variables. In addition to the total missing 
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transverse energy EJp^^^ and the transverse mass Mjjj used in the ATLAS analysis [36], we also 
include the often-used effective mass iJieS which is defined as the scalar sum of the four leading jet 
Pt's, the lepton px and E!^^^^. Signal and background events are generated using MadGraphS [37j, and 
showered in PYTHIA [38]. We use PGS |39j to perform the fast detector simulation, after modifying 
the code to implement the anti-/cj jet-finding algorithm with the distance parameter R = 0.4 [ID]. We 
simulated the events at 7 TeV center of mass energy so that we can cross check our results with the 
ATLAS paper [36j. |f| The signal production cross section is normalized to be the value calculated 
at NLO+NLL [41j, and the background production cross section for ti is normalized to be the value 
(Ttt{rn = 173 GeV, 7 TeV) = 163^519 pb, calculated approximately at NNLO [12]. In our studies, the 
leptonic decays of the top quarks contain r leptons. We adopt the same basic selection cuts on the 
objects in the final state as in Ref. [36] by requiring exactly one isolated electron or muon. 
^We focus on the light mass eigenstates, and calculate reach based on 100% decay to i + x'^. 

2n 



The transverse mass is defined by the formula Mt = u 2p!j,E^^^^[l — cos{(p^ — (p t )], where pj, is the pt of the 

leptons and 4>^ and (/f>^^ are the azimuthal angles of the lepton and E^^^^- 

■^We simulated 6 million events for the dileptonic and semileptonic ti backgrounds each. Although we only used the 
unmatched samples in this paper, we checked that a parton shower plus matrix element matched sample provided good 
agreement for the basic variable distributions with sufficiently large cuts. 
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Figure 1: The signal and background event distributions in three basic variables: E^^^^, Mt, and meg. 
The signal is 7 TeV production of a 500 GeV stop pair, each decaying to a top quark and a 100 GeV 
neutrahno. Ah the events in the plots have E!p^^^ > 150 GeV and Mt > 100 GeV. 



Fig. [T] shows the signal and background distributions in these three basic variables. For the signal 
we choose a stop mass of 500 GeV and the neutrahno mass of 100 GeV. We have included both 
dileptonic and semileptonic ti backgrounds. As one can see from the Mt distributions, the semi- 
leptonic tt background events mainly populate in the region with Mt < 150 GeV. Imposing a cut 
with Mt > 150 GeV will be an efficient way to suppress this background. We have also simulated the 
VF+jets background and found a similar distribution as the semi-leptonic ti background. With the 
Mt > 150 GeV cut, there is only a negligible number of the VF+jets background events left, so we 
will not include this background in what follows. The Mt cut is not effective at separating the signal 
and the dileptonic ti background events. On the other hand, cuts on E^^^^ and tjIcs can be used to 
significantly reduce this background, though it remains the biggest contamination in the direct stop 
production search. 

The diagram for the dileptonic ti background event topology is shown in Fig. [2l with dashed 
lines representing missing particles. Large E^™'^^ can arise due to the two missing neutrinos and the 
missing lepton. Also, the transverse mass Mt is not constrained by the W boson mass because of the 
additional missing particles. Because there are missing energies on both decay chains, the stransverse 
mass Mt2 [43^ |44] can be a natural variable to identify this type of background event. {Mt2 has 
been proposed to reduce ti and W~^W~ backgrounds in the di-lepton search channel [l5l 36].) The 
Mt2 for a given event can be interpreted as the minimal mother particle mass compatible with the 
postulated event topology and an assumed daughter particle mass [47]. The Mt2 is bounded from 
above by the mass of the mother particles in the decay chains if the assumed mass for the daughter 
particles is equal to (or less than) their true mass. By looking at the diagram in Fig. [2l we can 
define Mt2 and its generalizations or variations with the top quark as the mother particle for our 




Figure 2: The Feynman diagram for the ti background in the (dominant) dileptontic channel. The 
dashed lines represent missing particles at colhders, including a lost lepton that would otherwise 
exclude it as a background to our semileptonic stop signal. 



backgrounds. Our observables for the leading leptonic background are the 2 6-jets + one lepton 
+ E"™'*^*^ subsystem. In fact, the next-to-leading dominant semileptonic ti background also contains 
exactly the same subsystem if one disregards the jets from the W decay, so they may be used to bound 
this background too. On the other hand, the tt* signal has the additional missing energy source from 
the missing x particles. Consequently the corresponding variables can take larger values. 

In all Mr2-type variables, a minimization is performed over all possible ways of dividing E^^^^ 
between the two decay chains. More explicitly, the minimization is over all possible pairs of 4-momenta, 
each with an assumed mass, whose vector sum has transverse components that match E^^^^. The 
difference between variables comes in the assignment of visible and missing momentum to the two 
decay chains, along with invariant mass or Mj- constraints imposed on the hidden 4-momenta. In the 
following, we define three Mj'2-type variables with background endpoints roughly at the top mass. 
These new variables are not expected to be completely independent, so their performances will be 
evaluated in the next section. 

The first variable is basically the Mt2 of the tt — )• bW~^bW~ subsystem, which is denoted as M^g- 
Interpreted in the original Mj'2 context, it assumes a "missing on-shell W" on each side of the decay 
chain. Since the lepton momentum results from the W decay, we add it to the E^^^^. It is defined as 
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where the W mass is assigned for both pf and P2 and jet masses of pb^ and pb2 are calculated from 
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Figure 3: Schematic of M|.2, along with its signal and background event distributions. In standard 
Mt2 calculations, the Mt of each decay chain is calculated using a visible and invisible momentum. 
These are indicated by solid and dashed circles in the left panel. The "invisible on-shell VK" momenta, 
pi and p2) are scanned in the minimization calculation. The sum of the transverse components is 
taken to equal E^^^^ plus the visible lepton momentum's transverse component. All the events in the 
plot have £;??^"" > 150 GeV and Mt > 100 GeV. 



their four-vectors. A diagram illustrating this, along with signal and background distributions of M^2 
are shown in Fig. [3l By using the true W boson mass, M^2 is bounded between the top mass and the 
W gauge boson mass for the dileptonic tt background where this topology is appropriate. (The exact 
bound holds only for a perfect detector.) 

To select the two candidate 6-jets, we divide all events into three categories. The first category 
contains exactly two 6-tagged jets in the four leading jets of pT, and we can use Eq. ([1]) directly. For 
the second category containing exactly one 6-tagged jet, we choose the two leading non-6-tagged jets 
as the other b-jet candidate and take the smaller of the two M!^2^s. For the third category with zero, 
three, or four 6-tagged jets, we assume that the two candidate 6-jets are contained in the leading three 
jets and we ignore 6-tagging information. There are three different combinations, among which we 
take the smallest as the final value of M^2 ■ 

For the second variable, we do not add the observed lepton momentum to EJ^^^^. Instead we define 
an asymmetric M'r2 [ 481 SS] by combining the 4-momenta of the lepton and a b-jet into one effective 
particle. The missing neutrino on that side is treated as massless. On the other side, the visible 
particle is the other b-jet (with its mass calculated from its four- vector) , and the invisible particle is 



an on-shell W. This variable is defined as 
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The two ^-jet candidates are chosen by the same procedure as in the previous case. There are two 
ways to pair the lepton with one of the two 6-jets, and the combination which produces a smaller M!^2 
is chosen. A diagram illustrating the calculation, along with signal and background distributions are 
shown in Fig. HI 
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Figure 4: Schematic of Mpg, along with its signal and background event distributions. As compared 
to our previous variable M|,2, the 4-momentum of the visible lepton and 6-jet are combined together 
as one effective visible particle, and p2 is treated as an "invisible on-shell VF" when calculating Mt 
for that side. Again, all events in the plot have E^''^^ > 150 GeV and Mt > 100 GeV. 



Each of the two Mt2 variables defined above did not fully utilize the information available for the 
background event topology: the two intermediate W bosons are on-shell and one of them produces 
the observed lepton together with a neutrino. We can define a new kinematic variable as the minimal 
mother particle mass (the top quark mass in this case) which can be compatible with all the transverse 
momentum and mass-shell constraints of that topology for a given event. Here, the top quark mass is 
not explicitly used, only implicitly bounded by the event. This is in the same spirit as interpreting Mt2 
as the minimal mother particle mass compatible with the minimal kinematic constraints [47] except 
that all mass-shell constraints on the cascade decay chain are used|l|. One might expect to get a variable 



The mass-shell constraints are not sufficient to fully reconstruct each event. 



which is more sensitive to this background topology because of the additional kinematic information 
applied in the definition. Specifically, the variable M^^ (where the superscript W represents the on- 
shell intermediate W information is included when combining lepton and neutrino) can no longer be 
cast into the "maximum of two side's Mt" form, but is instead defined directly as the minimization |f| 
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Figure 5: Schematic of Mj^, along with its signal and background event distributions. Here all of the 
information is used, including the W^-on-shell mass condition on both sides. As with the other variables, 
P2 is the entire missing on-shell W , but pi is the neutrino that gets paired with the visible lepton to 
form the other on-sheU W. All the events in the plot have E^^'''' > 150 GeV and Mt > 100 GeV. The 
events with no compatible top mass under 500 GeV are placed in the last bin. 



The diagram, along with signal and background distributions are shown in Fig. O We use the same 
method as before to pick the two 6-jets, and a method similar to that for M^2 is used to choose which 
6-jet gets paired with the visible lepton. Calculating this variable can be done efficiently in a similar 
way as the Mt2 calculation in Ref. [47j by generalizing the method there to this case. For perfect 
measurements, this variable for the dileptonic ti backgrounds is less than the true top quark mass 
since the top mass should be compatible with all background events. On the other hand, the signal 
events do not need to satisfy such a bound, because of its different topology and additional missing 
massive particles x- For some of the signal events we may not even be able to find a compatible 
mass because we apply the variable to a wrong topology with the wrong mass-shell conditions. The 



The programs for calculating all new variables defined in this paper can be downloaded at 
https : //sites .google . com/a/ucdavis . edu/mass/ 



background distributions indeed lie mostly below the top quark mass, while a significant number of 
signal events have no solution below 500 GeV and they are included in the last bin. 

One can see from the plots in Figs. [3l [H and [5] that a cut on these variables around the top quark 
mass could be an effective way to suppress the main background. It is not clear a priori which one 
will have the best performance when the experimental smearing and detector resolution effects are 
taken into account, and whether there is still enough independent information among them so that a 
combination of them can give some further improvement. In the next section we will make a critical 
comparison of the performances of these variables and their combinations. 

3 Performances of New Kinematic Variables 

To quantify the power of these kinematic variables, we optimized a simple cut-and-count experiment 
involving three stop masses (400, 500, 600) GeV, with the neutralino mass being fixed at 100 GeV and 
100% branching ratio of stop decaying to top plus neutralino. We simulated the signal and background 
events at 7 TeV to compare with the existing ATLAS study. Although the LHC will run at 8 TeV 
this year, the relative performance of each kinematic variable will not be affected much. We will 
comment on the 8 TeV case in the next section. For each possible set of cuts, we use the NLO cross 
sections of stop pair productions multiplied by cut efficiencies to estimate the number of signal (S) 
and background (B) events expected in 20 fb^^ of data. The Poisson probability that pure background 
would fluctuate up to at least S + B events is given by ROOT's two-parameter F function, which can 
be evaluated at non- integer parameters: 

p = ^ -rre~^ = TMath :: Gainma(5' + 5,5). (4) 

k=S+B 

We translate this probability into a gaussian-equivalent significance (cr) in terms of standard devia- 
tions|£| This approaches S/yB for large signal and background, but by handling the small-number 
statistics, we avoid extreme cuts. By finding cuts that maximize this significance reach, we can 
estimate the power of of any set of kinematic variables. 

To evaluate the performances of the new kinematic variables defined in the previous section, we 
include them with a basic set of cuts on (E^^^^, Mt, w-off)- Within the basic set of variables, the 
E^^^^ is most powerful in discriminating the signal and the background. A cut on AIt > 150 GeV 
is imposed to remove the VF+jets background, and after this cut, the semileptonic ti background is 
virtually eliminated. We found that further increasing the Mt cut will hurt the signal significance. 



'This translation is also handled by ROOT as cr = TMath :: NormQuaiitile(l — p). 



Minimum Cuts 



m. 



stop 



400 GeV 



Eimiss 


rrics 


Mjf^ 


Mif, 


Mt2 


S20 fb-^ 


B20 fb-^ 


S/B 


a 


(150) 


- 


- 


- 


- 


129.3 


738.4 


0.17 


4.62 


202 


- 


- 


- 


- 


82.1 


208.8 


0.39 


5.34 


202 


491 


- 


- 


- 


81.7 


202.7 


0.40 


5.38 


202 


502 


- 


100 


- 


74.4 


147.9 


0.50 


5.66 


202 


502 


- 


- 


157 


55.9 


66.7 


0.84 


6.09 


200 


562 


177 


- 


- 


50.7 


48.5 


1.05 


6.33 


200 


564 


176 


95 


99 


50.0 


46.1 


1.08 


6.38 



Table 1: Cuts optimized for significance to discover 400 GeV stop and 100 GeV neutralino for 20fb~^ 
at 7TeV. Runs optimized over different cut variables are sorted by increasing significance (roughly 
S/vB for large numbers of events, where the improvement factor becomes independent of integrated 
luminosity.) All optimizations began with E^^^^ > 150 GeV (first row) and include a fixed Mt > 
150 GeV cut to eliminate VF+jets. The starting cuts yield 12120 simulated ti dileptonic and 24 
semileptonic background events, and 972 simulated signal events. The numbers of the signal and 
background events in the table are rescaled to 20 fb~^ luminosity. The Mt cut was for the Vl^+jets 
background, and increasing it beyond 150 GeV never helped with the ti backgrounds, so its column is 
not shown. 



so we fixed the Mt cut to be at 150 GeV in our study. The effective mass meg is not as useful as 
EJ^^^^, and its inclusion will only slightly improve the results. However, it is a simple variable and 
is widely used, so we still include it in the optimization. We added the new kinematic variables 
(M^2, M!^2j ^t2) °^^ ^^ ^ time to the basic set and compare improvements by re-optimizing the cuts 
on all variables. We also combined all three new variables with the basic set to see if there is any 
independent information among these new variables which can further improve the significance. The 
results are listed in Tables [l] and [2l 

Prom these results we see that indeed the new kinematic variables can improve the signal signifi- 
cance on top of the basic variables. (The effective mass rriefj, on the other hand, does not help very 
much.) For the 400 GeV stop, Mj^ has the best performance and M^g is the least useful one among 
the three variables. This is in accordance with our expectation of these three, as the variable M^ 
contains the most kinematic information of the background event topology while M^g contains the 
least. For a heavier stop of 500 or 600 GeV, the performances of the three variables are actually com- 
parable. (One should not take the small differences seriously due to the limited statistics.) These new 
variables are highly correlated and not much improvement can be gained by combining all of them. 
Other than improvement on the discovery sensitivity via S/\/B, we also note that from Tables [J and 
[2] the improvement on S/B is even more dramatic after including the new variables defined in this 
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Minimum Cuts 



m. 



stop 



500 GeV 



Eimiss 


rrics 


Mjf^ 


Mif, 


Mt2 


•520/6-1 


^20 fb-^ 


S/B 


a 


(150) 


- 


- 


- 


- 


34.0 


738.4 


0.05 


1.23 


303 


- 


- 


- 


- 


11.4 


16.6 


0.69 


2.49 


303 


659 


- 


- 


- 


11.4 


16.1 


0.70 


2.50 


299 


709 


172 


- 


- 


9.8 


6.2 


1.59 


3.19 


291 


743 


- 


163 


- 


7.9 


3.6 


2.21 


3.20 


300 


708 


- 


- 


170 


9.4 


5.6 


1.69 


3.20 


291 


742 
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123 


109 


9.0 


4.4 


2.04 


3.34 



Minimum Cuts 



'^stop 



600 GeV 



T7>miSS 
H/rp 


"leff 


M^ 


M|2 


Mt2 


•520/6-1 


^20 fb-^ 


S/B 


(T 


(150) 


- 


- 


- 


- 


16.7 


738.4 


0.02 


0.60 


377 


- 


- 


- 


- 


4.5 


3.0 


1.49 


2.04 


345 


696 


- 


- 


- 


6.1 


6.3 


0.97 


2.05 


337 


727 


168 


- 


- 


5.9 


3.0 


2.01 


2.66 


337 


726 


- 


- 


168 


5.8 


2.7 


2.17 


2.69 


333 


740 


- 


157 


- 


5.3 


2.1 


2.59 


2.73 


332 


741 


168 


148 


91 


5.5 


2.1 


2.67 


2.81 



Table 2: Cuts optimized for significance to discover 500 GeV and 600 GeV stops with 100 GeV neu- 
tralinos for 20fb~^ at 7TeV. Again, all runs began with EJ^^^^ > 150 GeV and include a fixed 
Mt > 150 GeV cut (not shown), where there are 2115 and 1938 simulated events for 500 GeV and 
600 GeV stops and the same number of background events as in Table [TJ Cuts on EJp^^^ and M^ still 
do almost as well as optimization over all variables, but here these additional cuts can improve S/B. 



paper. So, the systematic errors for the actual experimental searches can further reduced. 

We also tried a few small variations of these new variables and did not obtain better results. 
For example, in M^2! using zero mass for W, or not adding the lepton momentum to E^^^^ yields 
very similar results, and these variations are more than 95% correlated. Assuming the transverse 
momentum of the missing neutrino from the W decay is in the same direction as that of the observed 
lepton gives a worse result, since the W bosons in the background events are in general not highly 
boosted. An Mt2 variable motivated by the signal topology by combining one 6-jet with the lepton 
and the other 6-jet with two non-6-jets also does not help. Our results indicate that if one wants to 
choose a minimal set of variables for the semileptonic channel search of the stop direct production for 
a wide range of the stop mass, the set {EJfl^^'^, Mt, M^) (even without m^g) can achieve nearly the 
maximal discriminant power of combining many different variables. 
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4 Conclusions 

The LHC will be running at 8 TeV in 2012 and is anticipated to achieve a larger integrated luminosity. 
We expect an even higher mass reach for the stop search compared with the numbers obtained in the 
previous section based on 7 TeV. To estimate the exclusion or discovery sensitivity of the stop with 
a 20fb~^ luminosity, we calculated the stop signal cross sections at tree level using MadGraphS and 
applying the same i^-factor at the 7 TeV LHC to take into account the QCD NLO corrections. The 
same procedure is applied to the tt backgrounds to obtain the approximate NLO production cross 
section at the 8 TeV LHC. The total ti cross section is calculated to be 231.8 pb. With the help of 
the new kinematic variables discussed in this paper together with the basic variables EJ^^^^, Mt and 
rrieff, we found that for mi = 650 GeV and m^o = 100 GeV with Br(t — )■ t + x") = 100%, the stop can 
show up at the 4o" level if we ignore the systematic errors. The 95% C.L. exclusion reach can go up to 
around 700 GeV. If there is no excess found in the 8 TeV run, it will dent the hope of a non-fine-tuned 
SUSY solution to the hierarchy problem, unless the stop has some more exotic signatures in some 
non-standard scenarios, such as degenerate spectrum or i?-parity violation, etc. 

In this paper we have focused on suppressing the ti backgrounds for the search of direct stop 
production. However, the Standard Model ti production is a major background for a wide range 
of new physics searches at the LHC. The kinematic variables proposed here could also be useful in 
improving searches for other new physics where the ti constitutes the main background with a large 
missing transverse momentum for the signals. 

Comparing the performances of different variables shows that in general the more kinematic in- 
formation a variable contains, the more discriminant power it can possess. For more specific new 
physics searches where both the signal and main background event topologies are known, it is worth 
designing kinematic variables which carry as much information of the signal and/or the background 
events as possible to achieve the maximal discrimination between them. The strategy of constructing 
new kinematic variables discussed in this paper could be readily generalized to other cases. 
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